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Abstract We explore the implications of Bcrezinskii- 
Kosterlitz-Thoulcss (BKT) critical behavior and variable- 
range hopping on the two dimensional (2D) quantum 
superconductor-insulator (QSI) transition driven by tun- 
ing the gate voltage. To illustrate the potential and the 
implications of this scenario we analyze sheet resistance 
data of Parendo et al. taken on a gate voltage tuned 
ultrathin amorphous bismuth film. The finite size scal- 
ing analysis of the BKT-transition uncovers a limiting 
length L preventing the correlation length to diverge 
and to enter the critical regime deeply. Nevertheless 
the attained BKT critical regime reveals consistency 
with two parameter quantum scaling and an explicit 
quantum scaling function determined by the BKT cor- 
relation length. The two parameter scaling yields for 
the zero temperature critical exponents of the QSI- 
transition the estimates 3/2, z ~ 3 and V ~ 1/2, 

revealing that hyperscaling is violated and in contrast 
to finite temperature disorder is relevant at zero tem- 
perature. Furthermore, zv ~ 3/2 is also consistent with 
the two variable quantum scaling form associated with 
a variable-range hopping controlled insulating ground 
state. 
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1 Introduction 

Continuous quantum-phase transitions (QPT) are tran- 
sitions at absolute zero in which the ground state of a 
system is changed by varying a parameter of the Hamil- 
tonian [1,2 ,.3] . The transitions between superconduct- 
ing and insulating behavior in two-dimensional systems 
tuned by disorder, film thickness, magnetic field or with 
the electrostatic field effect are believed to be such tran- 
sitions [SEE]. 

Here we present a detailed analysis of the tempera- 
ture and gate voltage dependent sheet resistance data 
of Parendo et al. [51(7] . to explore in an ultrathin amor- 
phous bismuth film, the nature of the normal state to 
superconductor (NS) phase transition line including its 
end point at zero temperature, as well as the normal 
state to insulator (NI) crossover. At the common end- 
point of the two lines the film is supposed to undergo 
a quantum superconductor to insulator (QSI) transi- 
tion, separating the superconducting from the insulat- 
ing ground state. 

Considering the NS- transition line we explore the 
compatibility with gate voltage dependent Berezinskii- 
Kosterlitz-Thouless (BKT) critical behavior [8/9 ■ Our 
analysis of the temperature dependence of the sheet re- 
sistance at various fixed gate voltages uncovers a round- 
ed BKT-transition. The rounding is fully consistent with 
a standard finite-size effect whereupon the correlation 
length is prevented to grow beyond a limiting length 
L. Potential candidates for the limiting length include 
the failure in cooling [6,7, limited homogeneity due to 
local strain or a heat current. A nonzero heat current 
drives the system away from equilibrium and creates a 
temperature gradient which implies a space dependent 
temperature. Because the correlation length does not 
exhibit the usual and relatively slow algebraic diver- 
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gence as T c is approached, the BKT-transition is par- 
ticularly susceptible to such finite-size effects. On the 
other hand there is the Harris criterion [TOlHTj . stat- 
ing that short-range correlated and uncorrelated disor- 
der is irrelevant at the unperturbed critical point, pro- 
vided that v > 2/D, where D is the dimensionality of 
the system and v the critical exponent of the finite- 
temperature correlation length. With D = 2 and v = 
oo, appropriate for the BKT-transition [51(5], this disor- 
der should be irrelevant. In this context it is important 
to recognize that the existence of the BKT-transition 
(vortex-antivortex dissociation instability) in 4 He films 
is intimately connected with the fact that the interac- 
tion energy between vortex pairs depends logarithmic 
on the separation between them. As shown by Pearl [T^] 
vortex pairs in thin superconducting films (charged su- 
perfluid) have a logarithmic interaction energy out to 
the characteristic length X2D = A 2 /g?, beyond which the 
interaction energy falls off as 1 Jr. Here A is the magnetic 
penetration depth of the bulk. As A2D increases the 
diamagnetism of the superconductor becomes less im- 
portant and the vortices in a thin superconducting film 
become progressively like those in 4 He films [13] . Ac- 
cording to this \id >> min[W, L] is required, where 
W and L denote the width and the length of the per- 
fect sample. Invoking the Nelson-Kostelitz relation [2] 
\2D (T c ) = A 2 (T c ) Jd = <P 2 J (32ir 2 k B T c ) it is readily 
seen that for sufficiently low T c 's and min [W, L] « 1 
cm this condition is well satisfied. As a result a rounded 
transition uncovers a limiting length which is shorter 
than that resulting from the finite magnetic " screening 
length" A215. Nevertheless, for sufficiently large L the 
critical regime can be attained and a finite-size scaling 
analysis provides good approximations for the limit of 
fundamental interest, L —¥ 00 [15j . 

As will be shown below, the finite-size scaling anal- 
ysis uncovers rounded NS- phase transitions along the 
line T c (V g ), ending and vanishing at the critical gate 
voltage V gc where the QSI transition occurs. The crit- 
ical properties along this line are fully consistent with 
BKT-behavior subjected to a gate voltage dependent 
limiting length which decreases substantially by ap- 
proaching the QSI critical point at V gc . According to 
this and in agreement with previous studies [161117] . 
electrostatic tuning does not change the carrier density 
only but the inhomogeneity landscape as well. Further- 
more it is shown that the BKT-behavior leads to an 
explicit and extended expression of the standard quan- 
tum scaling form of the sheet resistance [ll[2l[3] 

The subscripts ± denote the branch resulting from the 
the Nl-crossover (+) and NS-transition (-), respectively 



c is a nonuniversal coefficient of proportionality, the 
gate voltage V g the tuning parameter with critical value 
V gc . z is the dynamic and V the critical exponent of 
the zero temperature correlation length, £ (T = 0, V g ) = 
£0 (T = 0) \5 — S c \~". G± (y) is a universal scaling func- 
tion of its argument such that G± (y = 0) = 1. The 
BKT-behavior yields with the variables Rq_ (V g ) and y 
the two parameter scaling form 



= G ± (y) . 



(2) 



R(T,V g ) 
tf0± (Vg) 

Noting that i?o- (V g ) exhibits a substantial gate volt- 
age dependence the one parameter scaling form ([1]) ap- 
plies asymptotically only (Rq- (V g ) —> R c )- As the film 
undergoes at finite temperature a BKT-transition the 
scaling function G_ (y) exhibits at the universal value 
y c a finite temperature singularity. The leading behav- 
ior of the phase transition line is then fixed by 



Tc(V a ) = -\V g -V gc \ 



(3) 



Contrariwise the BKT expression for the sheet resis- 
tance yields the explicit scaling function 



G_ (y) = exp (-b R (T/T c (V g ) - 1)- 1/2 ) 



(4) 



which applies for any T > T c . b R is a non- universal 
dimensionless constant. In this context it is essential 
to notice that the universal critical behavior is entirely 
classical for T c > [TS], including the characteristic 
form of the BKT correlation length which determines 
G_ (y). Quantum fluctuations enter via renormaliza- 
tionofi?o- {V g ), the constraint i? _ (V g ) — > i?o- (V gc ) — 
R c and T c (V g ). In contrast to the standard scaling form 
|T]) it applies for any T > T c (V g ) down to the critical 
endpoint of the BKT-line where Eq. ([3]) holds because 
T/T c approaches T/T c = y c /y. Using the finite size cor- 
rected estimates for i? - (V g ) and the BKT-line T c (V g ) 
we observe in an intermediate range of the scaling vari- 
able T/T c (V g ) a satisfactory collapse of the sheet re- 
sistance data onto the BKT scaling function (jj). The 
lower bound of this range stems from temperatures con- 
siderably above T c where BKT fluctuations no longer 
dominate and the upper bound is traced back to a fi- 
nite size effect preventing the correlation length to grow 
beyond a limiting length L(V g ). Taking the substan- 
tial gate voltage dependence of i?o- into account we 
observe close to the QSI-transition consistency with 
T c (V g ) oc \V g - V gc \ zV (Eq. ©) for zV ~ 3/2 and the 
corresponding expression for the BKT-scaling function 



(1) G_ (y) = exp (-b R (y c /y - 1) 1/2 ) 



(5) 



The estimate zv ~ 3/2 coincides with the value ob- 
tained from the analysis of the sheet resistance data 
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of gate voltage tuned epitaxial films of La2- 2; Sr :r Cu04 
that are one unit cell thick [171119] . 

To estimate the exponents z and 77 we note that ac- 
cording to the finite temperature dynamic scaling the- 
ory i?o_ (V g ) is proportional to ^ (V g ) where £o (V g ) 
is the amplitude of the BKT correlation length 

Z(T,V g ) = 

£o (V g ) ■ exp (- (b R /2) (T/T c (V g ) - iy 1/2 ) . (6) 

Neglecting the quantum constraint of a finite critical 
sheet resistance R c the amplitude i?o- (V g ) oc £g~ 2 (V^) 
should match the gate voltage dependence of the zero 
temperature counterpart £~ 2 (T = 0, V g ) which tends to 
zero at the critical gate voltage because £ (T = 0, V g ) di- 
verges as £ (T = 0, V g ) = Co (T = 0) (V g - V gc y v . The 
quantum constraint is then accounted for in terms of 

BSE] 

Ro-(V g ) - R c cx ^ 2 (V g ) oc C 2 (T = 0, V g ) 

<x(V a -V gc )- w <xJ^/'. (7) 

Using the estimates for Rq-(V 9 ), R c and T c (V g ) result- 
ing from the finite size scaling analysis of the temper- 
ature dependence of the sheet resistance at fixed gate 
voltages this relationship provides a direct way to mea- 
sure the critical exponents z and 77, and allows for a 
check on the value of zv. We obtain z ~ 3 and 77 ~ 1/2. 
Noting that D + z — 2 + z ~ 5 exceeds the upper criti- 
cal dimension D u = 4 the critical exponent of the zero 
temperature correlation length should indeed adopt its 
mean-field value V = 1/2. However the fate of this clean 
QSI critical point under the influence of disorder is con- 
trolled by the Harris criterion [TOlIll] . If the inequal- 
ity V > 2/D is fulfilled, the disorder does not affect 
the critical behavior. If the Harris criterion is violated 
(77 < 2/D), the generic result is a new critical point 
with conventional power law scaling but new exponents 
which fulfill 77 < 2/D and V ^ 1/2. Because disorder is 
relevant it drives the system from the mean-field to an 
other critical point with different critical exponents and 
with that to 77 ^ 1/2. Unfortunately the available data 
are too sparse to derive a more precise estimate of 77. 
Otherwise it is clear that D + z — 2 + z ~ 5 exceeds 
the upper critical dimension D u = 4 so that the the 
equivalence between quantum phase transitions in sys- 
tems with D spatial dimensions and the ones of classical 
phase transitions in (D + z) dimensions does not apply. 

Concerning the Nl-crossover it is known that in Bi- 
films spin-orbit coupling is important because its strength 
strongly depends on the atomic number Z and bismuth 
is the heaviest group- V semimetal [20] . In the presence 
of strong spin-orbit coupling, the spins rotate in the 
opposite direction and the interference is destructive, 
which leads a reduction of weak localization, known as 



weak anti- localization [2T] . Given the previous evidence 
for this behavior in Bi-films [20] one might expect that 
weak anti-localization controls the Nl-crossover. Other- 
wise it is also expected that in disordered films variable- 
range hopping controls the insulating phase in terms of 
a (T, V g ) = a h (V g ) exp (- (T (V g ) /T) 1/3 ) [22], trans- 
forming to the scaling function 

G+ (y) ~ exp(/y 1 / 3 ), (8) 

with the crossover temperature To (V g ) oc f 3 | V g — V gc \ zv 
In this context it should be kept in mind that the ir- 
relevance of disorder at the BKT-transition does not 
exclude its relevance at zero temperature. Indeed our 
estimate 77 ~ 1/2 violates the Harris criterion [TUlHTj 
and with that is the film considered here at T = sen- 
sitive to disorder. For this reason we assume that the 
insulating ground state is controlled by variable range 
hopping. We observe that the associated quantum scal- 
ing function (Eq. ©), combined with a strongly gate 
voltage dependent crh(V g ) — l/Ro+(V g ) mimics the 
sheet resistance data remarkably well, suggesting that 
the insulating phase of the Bi film considered here is 
controlled by variable-range hopping. 

Essential conclusions of our analysis include: As 
Ro-(V g ) = Ro(V g ) and Ro + (V g ) exhibit a strong gate 
voltage dependence is two parameter scaling with the 
variables i?o± (V g ) and y — c \V g — V gc \ zv /T unavoid- 
able to attain the QSI-transition. The piecewise data 
collapse onto the respective scaling function G± (y) 
extends over two decades of the scaling variable. The 
BKT-line leads to an explicit expression of the scaling 
function G_ (y) and the evidence for a variable-range 
hopping controlled insulating phase determines G+ (y). 
For large values of the scaling argument the deviations 
from these scaling functions are traced back to a finite 
size effect with a limiting length which increases with 
reduced gate voltage. In the lower branch G_ (y) devi- 
ations occur for small arguments as well because BKT 
fluctuations no longer dominate. The scaling function 
of the upper and lower branch are very different but 
both exhibit at the QSI-transition an essential singular- 
ity. Our estimates for the critical exponents at the QSI 
critical point are consistent with the following proper- 
ties: zv ~ 3/2 emerging from the BKT-line yields to 
a satisfactory data collapse onto both branches of the 
quantum scaling function G± (y). The estimates z ~ 3 
and 77 ~ 1/2 emerge from the gate voltage dependence 
of i?o— (V g ) — R c - The Harris criterion [TUlfTT] implies 
however that 77 = 1/2 is excluded, disorder is relevant 
at zero temperature but irrelevant an finite tempera- 
ture. Furthermore it is shown that the application of 
one parameter scaling form([T]) leads to rather a differ- 
ent z77 value because the regime where Ro-(V g ) ~ R c 
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applies is not attained, because Ro-(V g ) is a nonuniver- 
sal function entering the BKT-scaling form of the sheet 
resistance. Nevertheless its value at the QSI- transition, 
Ro-(V gc ) = R c is expected to be universal and given 
R c = c R h/4e 2 ~ c R ■ 6.45 kft [23lf24lf25] . c R is a dimen- 
sionless constant. Provided that hyperscaling applies it 
depends only on the universality class of the QSI tran- 
sition [23]. In the present case does hyperscaling not 
apply, because Z? + z = 2 + z~5 exceeds the upper 
critical dimension D u = 4. In any case our estimate 
R c ~ 19.35 kn is of order h/e 2 ~ 25.8 kfl. 

The nature of the two-dimensional QSI-transition 
has been intensely debated [4"ll5ll2~i] . The scenarios can 
be grouped into two classes, fermionic and bosonic. In 
the fermionic case the reduction of T c and the magni- 
tude of the order parameter is attributed to a combina- 
tion of reduced density of states, enhanced Coulomb in- 
teraction and depairing due to an increase of the inelas- 
tic electron-electron scattering rate [2"6"1[2~7] . The bosonic 
approach assumes that the fermionic degrees of freedom 
can be integrated out, the mean square of the order pa- 
rameter does not vanish at T c , phase fluctuations domi- 
nate and the reduction of T c is attributable to quantum 
fluctuations and in disordered systems to randomness in 
addition 24,25,28]. In the Bi-film considered here the 
consistency with finite size limited BKT behavior un- 
covers clearly the bosonic scenario: The BKT-line sepa- 
rates the phases of uncondensed and condensed Cooper 
pairs, while at zero temperature the QSI critical point 
separates the phase of condensed and uncondensed in- 
sulating Cooper pairs. 

In Sec. II we sketch the theoretical background and 
present the detailed analysis of the resistivity data of 
Parendo et al. [HUT] . We close with a brief summary and 
some discussion. 



2 Theoretical background and data analysis 

2.1 BKT transition 

To explore the compatibility with BKT critical behav- 
ior we invoke the characteristic temperature dependence 
of the correlation length above T c [8,9 , 



Z(T,V g ) = Z (V g )exp 



b_R (V g ) 
2{T-T c {V g )f 2 1 



(9) 



Notice that for any T c > the universal critical behav- 
ior is entirely classical [15] . including the characteristic 
form of the BKT correlation length. Contrariwise £o ($) > 
bu (5) and T c (8) are subjected to quantum fluctuation 
renormalization and depend on the tuning parameter 8. 



bji is related to the vortex core energy and £o to the vor- 
tex core radius [16,29,30,31 . Invoking dynamic scaling 
the sheet resistance R scales in D = 2 as .2 



RcxC 



(10) 



where z c i is the dynamic critical exponent of the classi- 
cal dynamics. z c i is not questioned to be anything but 
the value that describes simple diffusion: z c i = 2 [3"!Z1 
133] . Combining these scaling forms we obtain 

2 



R(T,6) = / fr, {V g ) 
Ro-(S) \t(T,V g ) 



exp 



b R {V g 



(T-T c (V g )) 



1/2 



(11) 



The compatibility of sheet resistance data with the char- 
acteristic BKT behavior can be explored in terms of 

d\n(R(T,V g ))Y 2/3 



dT 



2/3 



(T-T c (V g )). 



(12) 



M(v g ) 

However, supposing that there is a limiting length L 
preventing the growth of the correlation length the tran- 
sition temperature T c cannot be approached and with 
that the attainable critical regime is reduced. The ex- 
tent of the attained BKT-regime can be extracted from 
a finite size scaling analysis. It implies that the mea- 
sured R (T, L) scales as [15] 

2 

= f{x), (13) 



R(T,L) 



£(T,L = oo) 



Z(T,L) 



R(T,L = oo) 
where 

x = Rr 1 (T, L = oo) oc (T, L = oo) 



(14) 



R (T, L — oo) and £ (T, L — oo) denote the respective 
behavior of the infinite and homogeneous system given 
by Eqs. ^ and (fTT| . respectively. f(x) is the finite size 
scaling function adopting the limiting behavior: 



1 :£(T,oc)«L 
g (L) x : £ (T, oo) »L 



(15) 



Accordingly, BKT-behavior, R (T, L) = R(T,L = oo) 
can be observed as long as £ (T, oo) < L , while for 
£ (T, oo) > L the scaling function approaches the finite 
size dominated regime where 

R{T,L) ong(L) oc L~ 2 . (16) 

Note that g (L) fixes the width of the intermediate re- 
gime where BKT-bchavior is observable. It shrinks with 
increasing g (L), i.e. with decreasing limiting length L. 
This power law neglects the possible occurrence of mul- 
tiplicative logarithmic correction i.e. R(L) oc g (L) oc 
L~' 2 1 In (L/&o)-[31] Noting that the finite size scaling 
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T(K) 

Fig. 1 R vs. T for various gate voltages V s for the 10.22 A 
film taken from Parendo et al. [5][7] for various gate voltages 

plot yields g(L), this correction leads to a saturation 
of the sheet resistance at low temperatures as well and 
that the data considered here is too sparse to uncover 
it, we neglect this correction. The occurrence of BKT- 
behavior also requires that T c is sufficiently below the 
mean- field transition temperature T c q. It can be esti- 
mated with the aid of the Aslamosov-Larkin (AL) ex- 
pression for the conductance [35] 

a = cr„ + cr / In (T/T c0 ) , (17) 
with 

2 

a = —r - 0-0152 kfT 1 , (18) 
8h 

where Gaussian fluctuations are taken into account. 

We are now prepared to explore the occurrence of 
BKT-behavior emerging from the sheet resistance data 
of Parendo et al. [Sl[7] depicted in Fig. [T] for an amor- 
phous 10.22 A thick Bi film taken at various gate volt- 
ages. It is apparent that most sheet resistance curves 
exhibit a rounded normal state to superconductor tran- 
sition. Nevertheless a reduction of the smeared tran- 
sition temperature with decreasing gate voltage can 
be anticipated. Accordingly BKT-behavior is expected 
in an intermediate temperature regime above T c only, 
and with that . To uncover this regime we proceed 
as follows: Invoking Eq. (fT2"|) consistency with BKT- 
behavior is established in terms of an intermediate tem- 
perature regime revealing a linear temperature depen- 
dence for the adjusted parameters (2/bn(V g )) 2 ^ 3 and 
T c {V g ). The amplitude Rq (V g ) is then estimated by ad- 
justing Eq. (1111) in the intermediate temperature regime 
to the measured temperature dependence of the sheet 
resistance with given (2/bn(V g )) 2 ^ 3 and T c (V g ). The 
quality of the resulting estimates for (2/bn (V g )) 2 ^ 3 , 
T c (V g ) and Rq (V g ) as well as the evidence for finite size 



limited BKT behavior is finally clarified in terms of the 
finite size scaling relation (IT31 . Results of this analysis 
are shown in Fig. [2^ uncovers consistency with BKT- 
behavior (HU), i.e. for V g = 19.5 V with (2/b R ) 2/3 ~ 
5.3 K" 1/3 and T c ~ 0.027 K in the intermediate regime 
extending from T ~ 0.055 K to T = 0.08 K. Above this 
regime BKT-fluctuations no longer dominate, while be- 
low its lower bound a pronounced finite size effect oc- 
curs. These features are well confirmed in Fig. [2b show- 
ing R vs. T. Indeed, the curves corresponding to Eq. 
(fTTj) fit the data in the intermediate temperature regime 
with the adjsuted R = 33.5 kfi, (2/b R ) 2/3 ~ 5.3 K" 1 / 3 
and T c ~ 0.027 K rather well. Furthermore the BKT- 
transition temperature is considerably below the mean- 
field counterpart as illustrated for V g = 19.5 V where 
T c ~ 0.027 K compared to T c0 ~ 0.04 K. Considering 
the finite size scaling plot in Fig. [5J; it becomes clear 
that the rounding of the transition is fully consistent 
with a standard finite size effect, preventing the corre- 
lation length to grow beyond a limiting length L. Not- 
ing that the intermediate regime where BKT-behavior 
is observable shrinks with increasing g, the rather large 
value .g = 9 kfi at V g = 19.5 V is consistent with the 
severely rounded transition. Nevertheless, in agreement 
with Fig. [2^l and Fig. [2Jd there is a window left where 
BKT-behavior occurs. 

The estimates for the BKT parameters and their 
gate voltage dependence are shown in Fig. [3] Fig. [3^ 
uncovers a BKT-transition line T c (V g ) with a quan- 
tum critical cndpoint around V gc — 10.5 V. Though 
these transition temperatures are not attained, the gate 
voltage dependent limiting lengths turned out to be 
sufficiently large to estimate them and the other pa- 
rameters of interest from the intermediate temperature 
regime where consistency with BKT-behavior was es- 
tablished and confirmed in terms of a finite scaling anal- 
ysis. The gate voltage dependence of g (L) reveals that 
electrostatic tuning does not change the carrier density 
only but the inhomogeneity landscape as well. Indeed 
g (L) oc L~ 2 increases from V g = 33 V to V g = 12.5 V 
by almost an order of magnitude. Note that g (L) fixes 
the width of the intermediate regime which shrinks with 
increasing g (L), i.e. with decreasing limiting length L. 

The solid and dashed lines in Fig. [3^ are 

T c (Vg) = 0.0011 KV- 3/2 • (V g - V gc ) 3/2 , (19) 
and 

T c (V g ) = 0.0065 KV- 2/3 (V g - V gc ) 2/3 , (20) 

with V gc = 10.5 V. They indicate limiting and effective 
power law behavior of the BKT phase transition line 
near its quantum critical endpoint where T c vanishes. 
In Fig. [5b we show the estimates for the T c dependence 
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a: 



1 / R(T,L=oo) (kn" 1 ) 

Fig. 2 a) (dln(R)/dT)- 2/3 vs. T derived from the data of 
Parendo et al. [S][T]. The lines are fits to Eq. I|12l) yielding esti- 
mates for T c and (2/b R ) 2 / 3 . (T c = 0.027 K and (2/&r) 2/3 = 
5.3 K" 1 / 3 .) b) R vs. T for the same gate voltages. The lines 
are Eq. © with T c , b R derived from Fig. 2a and adjusted _Ro 
to match the respective sheet reistance data above T c . The 
solid line is a fit to the Aslamosov-Larkin expression (|f 711 
yielding cr„ = 0.044 kfi -1 and T c0 = 0.04 K compared to the 
corresponding BKT transition temperature T c = 0.027 K. c) 
Finite size scaling plot R (T, L) /R(T, L = oo) vs. 1/R{T, L = 
oo). The solid line marks the occurrence of BKT-behavior 
and the dash-dot line the finite size dominated regime where 
R (T, L) ocg(L) oc L~ 2 . 



10 15 20 25 30 35 40 45 
V (V) 

Fig. 3 Gate voltage and T c dependence of the BKT- and 
finite size scaling parameters T c , Rq, bu and g/L 2 . a) T c and 
g vs. V g . The solid line is T c = O.OOlf KV" 3 / 2 (V g - V gc ) 3/2 
and the dashed one T c = 0.0065 KV" 2 / 3 (V g - V gc ) 2/3 with 
V gc = 10.5 V. b) Rq and b R vs. T c . The solid line is R = 
Rq c + 128T, 2/3 kQ with i? 0c = 19.35 kQ and the dashed one 
b R = 0.99T C 1/2 K 1 / 2 . c) _R vs. V g . The solid line is Rq = 
R c + lA{V g - V gc )kQ with V gc = 10.5 V and R c = 19.35 
kfi. 
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of i?o- and bji. i?o- is the amplitude of the BKT ex- 
pression for the sheet resistance (Eq. (fTTj) 1 ) and b R is 
related to the vortex core energy and should scale as 
fl6ll29l30l 



b R {V g )=b R TV\V g ) (21) 

A glance at Fig. reveals that this behavior is very 
well confirmed with the nonuniversal parameter 

b R ~ 0.99. (22) 

The estimates for i?o- (V g ) shown in Figs. [3Jd and c 
also reveal that this amplitude exhibits a pronounced 
T c and gate voltage dependence. In the limits T c — » 
and V g — > V gc ~ 10.5 V it tends to the critical value 
R 0c ~ 19.35 kfi. The dashed line in Figs. [3b is 

i?o- (T c ) = i? c + 128 kQK- 2/3 ■ T c 2/3 , (23) 

and the solid line in Figj3]: 



i? - (Vg) =R C + 1-35 kfW- 1 • (Vg - Vg, 



(24) 



They indicate power law behavior near the quantum 
critical point at V gc ~ 10.5 V and R c ~ 19.35 Ml where 
T c vanishes. The reliability of these relationships can be 
checked by eliminating i?o- (Vg) — R c and solving for 
T c (V g ). This yields Eq. (J2J) as it should be. 

An essential implication of relation ([2"Tj) is that the 
BKT expression (jTTJ) for the sheet resistance adopts the 
two variable scaling form 



R(T,Vg) = / CO (Vg) 

i?o-(K,) U(T,^) 



6 fl 



exp I — . (25) 

V (T/T c (V g )-l) 1/2 J 

Indeed it depends on Rq_ (V g ) and T c (V g ) only. Ac- 
cordingly, given R (V g ) and T c (V g ) the sheet resis- 
tance data R (T, V g ) plotted as R (T, V g ) /R - (V g ) vs. 
T c (V g ) JT should fall on a single curve given by the 
right hand side of Eq. (|2"5"j) . Deviations from this curve 
are expected close to T c (V g ) /T < 1 due to the finite 
size effect and for T c (V g ) /T « 1 where BKT fluctua- 
tions no longer dominate. This scaling plot is shown in 
Fig. 01 The piecewise collapse of the data onto the scal- 
ing function, marked by the solid line clearly confirms 
that the BKT regime is attained. For large values of the 
scaling argument it is bounded by the finite size effect 
with a limiting length which increases with reduced gate 
voltage (Fig. [3^) and for small arguments they reflect 
the fact BKT fluctuations no longer dominate. 

Evidence for a gate voltage driven QPT emerges 
from the phase transition line T c (V g ) pointing to an 
endpoint where T c vanishes (Fig. |3^) at a critical gate 
voltage Vg C ~ 10.5 V and the amplitude Rq (V g ) of 
the sheet resistance tends to the critical value Rq c ~ 




Fig. 4 Scaling plot R(T,V a )/R -(V g ) vs. T c (V g )/T for the 
data shown in Fig. fusing the T c (V g ) and Ro-(V g ) values 
shown in Fig. [3^ and Fig. [3J;, respectively. The solid line is 

exp {-b R / {T/T c (V g ) - 1) 1/2 ) in to Eq. ([25) with b R ~ 0.99 

(Eq. p2]l). 



19.5 k£l (Fig. [3fc). A characteristic property of a QPT 
at the endpoint of a BKT phase transition line follows 
from Eq. ((25]) by considering the isotherms. Indeed, as 
T c (V g ) tends to zero by approaching V gc at fixed T, the 
amplitude of the sheet resistance R(T, V g ) tends to the 
temperature independent value R c at V gc . Accordingly, 
the isotherms merge at the endpoint of a BKT tran- 
sition line T c (V g ) which is a characteristic property of 
a QPT. In Fig. [5j showing sheet resistance isotherms 
derived from the data depicted in Fig. [TJ we observe 
that the flow of the measured data points ( O? to a 
merging point near V gc ~ 10.5 V and R c ~ 19.35 kSl is 
well confirmed and consistent with respective values (*) 
derived from Eq. (1231) and the T c (V g ) and Ro(V g ) values 
shown in Fig. and Fig.[3t. In contrast, the solid and 
dashed lines, obtained from Eq. (|25p and the power law 
expressions for T c (V g ) and R (V g ) (Eqs. (JTSJ) and l[2"3|)) 
apply closer to the QPT only. 



2.2 Quantum phase transition 

The scaling theory of QPT yields for the sheet resis- 
tance to the one variable scaling form [1 , 23 



R(T,Vg 

Rr 



G± (y),y 



c\V g -Vgc\' 
T 



(26) 



c is a nonuniversal coefficient of proportionality, R c the 
critical sheet resistance, z the dynamic and 77 the criti- 
cal exponent of the zero temperature correlation length, 



£(T = 0,V/ 9 )=£ (T = 0) \V g - V gc \~ 



(27) 



G± (y) is a universal scaling function of its argument 
such that G± (y = 0) = 1. The subscript (+) marks 
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Fig. 5 Sheet resistance isotherms R{V g ) at T = 0.075 K 
and T = 0.05 K. O an d Q denote the measured data points 
derived from Fig. [T] * mark the values derived from Eq. I|25l) 
using the T c (V g ) and Ro(V g ) values shown in Fig. [3^i and 
Fig- [St- The solid and dashed lines are derived from Eq. I|25l) 
using the T c (V g ) and Ro(V g ) given by the power laws H19I) 
and 



the Nl-crossover and (— ) the NS-transition line. Due to 
this line, corresponding in the present case to the BKT- 
transition line, the scaling function G_ (y) exhibits at 
finite temperature a singularity at the universal vale 
y c of the scaling variable. The leading behavior of the 
phase transition line dill) is then fixed by 



T c (V g ) = -\V g -V gc r 

Noting that in this limit 

V _ T c (V g ) 



T 



(28) 



(29) 



the quantum scaling function G_ (y) adopts according 
to Eq. ([25]) the explicit BKT form 

G^(y) = cxp(-b R (y c /y-l)- 1/2 y (30) 

Because R c is fixed there is one scaling parameter only, 
namely y oc T c (V^) jT. Taking the gate voltage depen- 
dence of i?o into account we obtain the two parameter 
scaling form 

R(T,V g ) 



G- (y). 



(31) 



Ro- (v g ) 

The occurrence of a QPT at the endpoint of the 
BKT-transition line implies that the isotherms merge 
at the critical value of the tuning parameter V gc with 
sheet resistance R c . where at T = the QSI-transition 
occurs. In the BKT scenario where the NS-branch is 
considered this crossing point corresponds to the merg- 
ing point seen in Fig. [5] In the quantum scaling ap- 
proach the scaled data R (T, V g ) /R c vs. \V g - V gc f jT 
should then fall onto two branches crossing at the QSI 



critical point. The upper branch stems from the NI- 
crossover and the lower one from the NS- transition, zv 
is usually determined by minimizing the error in data 
collapse [4,5,6j[7]. The dashed line in Fig. [3^,, given by 
Eq. (l20l) . suggests that the effective exponent zv ~ 2/3 
mimics the overall gate voltage dependence of T c rea- 
sonably well. The resulting scaling plot, including the 
data shown in Fig.Q]is depicted in Fig. [5^,. Admittedly 
the data collapse is rather poor. For large values of the 
scaling argument the collapse is affected by the finite 
size effect with a limiting length which increases with 
reduced gate voltage (Fig. [3^) and for small values of 
the scaling argument one enters T » T c (V g ) where 
BKT fluctuations no longer dominate. More alarming 
are the large deviations from the BKT-behavior indi- 
cated by the solid line. This discrepancy is attributable 
to the fact that at this gate voltages the amplitude Rq- 
decreases substantially and differs considerably from 
the critical value R c . On the other hand adopting the 
suggestion of Parendo et al. the quality of the data 
collapse can be improved by restricting the data, i. e. to 
the interval 0.06 K < T < 0.1 K. As shown in Fig. [6b 
this improves the quality of the collapse considerably 
and suggests consistency with zv ~ 2/3. Although on a 
first glance compelling, this procedure is misleading be- 
cause the asymptotic regime where the quantum scaling 
form ([21))) applies is not attained as the gate voltage de- 
pendence of T c (see Fig. [3^) and the amplitude i?o- (V g ) 
reveals (Fig. [3b). 

Indeed, Fig. [3feo shows that z~P ~ 2/3 is at best an 
effective exponent failing near criticality where zv cer- 
tainly exceeds one and zv ~ 3/2 is consistent with the 
available data. Furthermore zv ~ 3/2 is consistent with 
the limiting behavior of the sheet resistance isotherms 
close to the quantum critical point as depicted in Fig. [S] 
According to Fig. [7] this value also leads to a piecewise 
excellent data collapse onto the BKT scaling function 
([3"0)) using the two parameter scaling form (j3"Tj) . For 
large values of the scaling argument the deviations stem 
from the finite size effect with a limiting length which 
increases with reduced gate voltage (Fig. [3^,) and for 
small arguments they reflect the fact that BKT fluc- 
tuations no longer dominate. Considering Fig. [3^ we 
observe that the attainable quantum critical regime ex- 
tends up to V g ~ 19.5 V and according to Fig. [5b up to 
the temperature T ~ 0.07 K. A restriction of the sheet 
resistance data, i.e. to the interval 0.06 K < T < 0.1 K 
implies then that in the scaling plot shown in Fig. [BJd an 
essential part of the quantum critical regime is missing. 

Next we turn to the gate voltage and T c dependence 
oft the BKT amplitude Rq shown in Figs.^js andUJ:. In 
D = 2 the sheet resistance, while diverging in the in- 
sulating phase, and zero in the superconducting phase, 
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Fig. 6 Quantum scaling plot R(T,V g )/R c vs. 
(V g -V g c) z ~/T with zv = 2/3, T c = 19.5 kQ and 
V gc — 10.5 V. a) For the data shown in Fig.[T] The solid line 
is Eq. J25) with R (V g ) = R {V gc ) = R c = 19.35 kf2,& R = 
0.99, and T c (V g ) = 0.0065 KV" 2 / 3 (V g - V gc ) zv with 
zv = 2/3 (Eqs. (gb} and (J2SJ) indicates the BKT behavior, 
b) For the same data but restricted to 0.06 K< T < 0.1 K 
and the same solid line. 



a: 




0.0-1 , 1 , 1 . , . 1— 

200 400 600 800 

(V g -V gc ) 3 ' 2 / T (V m l K) 



Fig. 7 Scaling plot R (T, V g ) /R - (V g ) vs. {V g - V gc ) zv IT 
with zv = 3/2, T c = 19.5 kQ, V gc = 10.5 V and The solid 
line is the BKT scaling function Eq. pOl) with b R = 0.99, 
T c (V g ) = 0.0011 KV" 3 / 2 (Vg - V gc ) z17 and zv = 3/2 (Eqs. 
(IT9l and dgHJ). fl = R -(V g ) is taken from Fig.[3E 



has a finite constant value R c right at the quantum 
critical point. Indeed along the BKT phase transition 
line approach both Ro-(V g ) and Rq_(T c ) the critical 
value R c ~ 19.5 kO. We observe that In the quan- 
tum critical regime outlined above (V g < 19.5 V and 
T c ;$ 0.07 K) the approach to R c is consistent with 
the power laws and Note that in analogy 

to 6fl(V^) and T c (V g ) the amplitude i? -(Vg) is also 
renormalized by quantum fluctuations. To estimate the 
exponents z and V we note that according to classical 
dynamic scaling Rq_ (V g ) is proportional to ^ 2 (V g ) 
where (V g ) is the amplitude of the BKT correlation 
length 

i(T,V g ) = 

Ho (V g ) exp (- (b R /2) (T/T c (V g ) - \y l/2 ) . (32) 

At low temperatures where T c (V g ) oc (V g — V gc ) zv (Eq. 

) applies it has to match the quantum scaling coun- 
terpart 

H{T,V g )<x{V g -V gc )- v f-(y) 

y = cT/ (V g - V gc ) zV oc T/T c , (33) 

extending Eq. ([2T]) to finite temperatures. The scaling 
function /_ (y) adopts the limits: 

The matching yields an explicit form of the scaling func- 
tion /_ (x) and reveals that £o (V g ) scales as {Vg) — 
£ (T = 0, V g ) oc (V g - V gc )~ v oc T~ x l\ On the other 
hand there is the quantum constraint of a finite criti- 
cal sheet resistance R c at T = which is incompatible 
with the classical dynamic scaling prediction i? - {V g ) 
oc £q 2 (V g ) because ^ 2 (V^) tends to zero as (V g — V gc ) 1 
The quantum constraint is then accounted for in terms 
of QMI] 

Ro-(V g )~R c ^l/e (T = Q,V g ) 

<x(V g -V gc ) 217 <xTy*. (35) 

Using the estimates for Ro_(V g ), R c and T c (V g ) result- 
ing from the finite size scaling analysis of the temper- 
ature dependence of the sheet resistance at fixed gate 
voltages this relationship provides a direct way to mea- 
sure the critical exponents z and 17, and allows for a 
check on the value of zv. Previously it was confirmed 
in the analysis of the sheet resistance data of the su- 
perconducting LaAlOs/SrTiOa interface [TB] and thin 
La2- K Sr ;E Cu04 films [17] . Furthermore aT c was 
verified in submonolayer superfluid 4 He films |31) . Com- 
bining Eqs. (|35p . (|23p and (|24p we obtain for the critical 
exponents of the QSI transition the estimates 

zv~ 3/2, z ~ 3,17 ~ 1/2. (36) 
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Noting that D + z = 2 + z~5 exceeds the upper critical 
dimension D u — 4 the critical exponent of the zero tem- 
perature correlation length should adopt its mean-field 
value v — 1/2. However the fate of this clean QSI criti- 
cal point under the influence of disorder is controlled by 
the Harris criterion [TUllllj . If the inequality V > 2/D 
is fulfilled, the disorder does not affect the critical be- 
havior. If the Harris criterion is violated (p < 2/D), 
the generic result is a new critical point with conven- 
tional power law scaling but new exponents which ful- 
fill V < 2/D and z7 ^ 1/2. Because disorder is relevant 
it drives the system from the mean-field to an other 
critical point with different critical exponents and with 
that to v ^ 1/2. Unfortunately the available data are 
too sparse to derive a more precise estimate of V. Other- 
wise it is clear that £> + z = 2 + z~5 exceeds the upper 
critical dimension D u = 4 so that the the equivalence 
between quantum phase transitions in systems with D 
spatial dimensions and the ones of classical phase tran- 
sitions in {D + z) dimensions does not apply. 



2.3 Quantum superconductor to insulator transition 

In disordered 2D systems weak localization manifests 
itself due to the wave nature of electrons since inter- 
ference effects double the classical probability for self- 
crossing paths and it slightly suppresses the conduc- 
tance 36 . In the presence of strong spin-orbit coupling, 
the spins rotate in the opposite direction and the inter- 
ference is destructive, which leads to an increase in the 
conductance, this effect being opposite to the weak lo- 
calization and known as weak anti- localization [21j . The 
influence of spin-orbit coupling increases for heavy ele- 
ments since it strongly depends on the atomic number 
Z, thus bismuth, being the heaviest group-V semimetal, 
is a suitable material for the study of its electronic 
transport properties through the analysis of the weak 
anti- localization contribution |20) . The fermionic part 
comprises the sheet conductance and its modifications 
due to weak localization, electron-electron- and spin or- 
bit interaction, yielding [37 

c„ = <J\ + cr 2 In (T) . (37) 
Note that 
e 2 

o- 2 = — ~ 1.23 x 1CT 2 kfT 1 (38) 
Tin 

is generally attributed to the electron-electron interac- 
tion contribution [38] . In the vicinity of the normal state 
to superconductor transition there is also a bosonic con- 
tribution due to Cooper pair fluctuations. In the Gaus- 
sian approximation it is given by the Aslamosov-Larkin 



(AL) expression (|17[) . In Fig. [5] we show the sheet con- 
ductance a vs. ln(T) for gate voltages covering the nor- 
mal state to superconductor transition and the normal 
state to insulator crossover. Considering the behavior 
at V g = V and V g — 7 V we observe consistency with 
weak anti-localization. Indeed the gate voltage depen- 
dent values for 02 are substantially smaller than the 
electron-electron interaction contribution [38^ . The sat- 
uration above ln(T) « —3 is again attributable to the 
finite size effect preventing the diverging length associ- 
ated with weak localization to grow beyond the limit- 
ing length L. Although there is evidence for weak anti- 
localization at V g = 12.5 V and even at V g = 11 V below 
ln(T) w —2 Cooper pair fluctuations dominate at low 
temperatures as the comparison with the AL- behav- 
ior (Eq. (|T7|) ) and BKT behavior reveals (see FigJ5b). 
Here the approach to a normal state to superconductor 
transition, limited by the finite size outlined above, is 
apparent. Noting that at T c , a oc £ 2 oc L 2 , the gate volt- 
age dependence of the saturation in a reflects the reduc- 
tion of the limiting length L with reduced gate voltage 
(see FigJSJi). Similarly, as the gate voltage increases on 
the insulating side from V to 7 V the saturation level 
decreases because the diverging length is limited again 
by a length which decreases by approaching the criti- 
cal gate voltage V g ~ 10.5 V from below. The features 
emerging from Fig. [T] can be summarized as follows: 
A characteristic feature of the normal state appears 
to be the competition between weak anti-localization 
and superconductivity. For V g < 10.5 V superconduc- 
tivity is suppressed and a finite size limited insulat- 
ing ground state is approached. This approach is in 
agreement with previous work |20) consistent with weak 
anti-localization. On the contrary for V g > 10.5 V pair 
fluctuations dominate at sufficiently low temperatures 
and a finite size limited approach to the superconduct- 
ing ground state is observed. From the merging point 
seen in Fig. [5] we know that the flows to the supercon- 
ducting and insulating ground states are separated by 
the the quantum critical point at V gc — 10.5 V and 
cr c = l/i? c — 0.0517 kfl . The temperature depen- 
dence of the sheet conductance at various gate voltages 
seen in Fig. [5£i is then consistent with the approach to 
the QSI-critical point with finite critical sheet conduc- 
tivity. 

Otherwise it is also expected that in disordered films 
variable-range hopping controls the insulating phase 
whereupon the sheet conductivity scales as [22] 

a (T, V g ) = u h (V g ) exp (- (T„ (V g ) /T) 1/3 ) . (39) 

Tq (V g ) denotes the variable-range hopping crossover 
temperature. In this context it should be kept in mind 
that the irrelevance of disorder at the BKT-transition 
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Fig. 8 Sheet conductance a (T, V g ) vs. ln(T) for data shown 
in Fig. ^ The horizontal thin line marks the critical sheet 
conductance a c = l/R c ~ 0.0517 kQ _1 . The dashed line 
is Eq. (I37|) describing the fermionic contribution with o\ = 
0.062 kfi" 1 and cr 2 = 4.69 x 10" 3 kQ-^nd the solid one 
with cri = 0.063 kfi- 1 and cr 2 = 6.01 x 10" 3 kfi -1 . The 
dotted and dash-dot-dot lines are the AL expression (|17l) 
describing the gaussian bosonic contribution to the normal 
state to superconductor transition with <r n = 0.0514 (kfi -1 ), 
T c0 = 0.036 K, <r„ = 0.0505 kQ" 1 and T c0 = 0.028 K. The 
dash-dot line is the BKT behavior resulting from Eq. (|25l) 
with b R = 0.99 (Eq. R - = 22.05 kfi and T c = 0.002 K 

given in Fig. [3] 



does not exclude its relevance at zero temperature. As 
our estimate v ~ 1/2 violates in D = 2 the Harris crite- 
rion we know that disorder is relevant at zero temper- 
ature. For this reason we consider below an insulating 
ground state controlled by variable-range hopping. An 
approach to the QSI critical point requires then that 
<Jh (V g ) -> a h (V gc ) = a c and T (V g ) -> T (V gc ) = 0. 
To explore this scenario quantitatively we note that Eq. 
(j31?)) leads according to Eq. (|2l>| to the quantum scaling 
form 



flot «) » (V.) 
with the explicit scaling function 

G+ (y) ~ expt/y 1 / 3 ), 

and the crossover temperature 

w) = / 3 iK,-K, c r, 



T 



(41) 



(42) 



where the value of zv applies on both side of the QSI- 
transition. In Fig. [5^ we depicted the scaling plot 
R(V g )/R 0+ (V g ) vs. (T (V g )/T) 1/3 by choosing 
i?o+ (V g ) = \joh (V g ) and To (V g ) to minimize the error 
in data collapse onto the scaling function 
exp(T (V g ) /T) 1/3 . Though the quality of the collapse 
is satisfactory and points to a variable-range hopping 
controlled insulating ground state, it is again subjected 



to the limiting length identified with the finite size scal- 
ing analysis of the BKT transition. It leads to the sat- 
uration at low temperature. Additional insight emerges 
from the gate voltage dependence of the adjusted pa- 
rameters To (V g ) and i?o+ (V^) shown in Figs. |9t> and 
[5J:. In contrast to the BKT-transition line (see Fig. [3^) 
the quantum scaling form (I42[) of the crossover temper- 
ature To (V g ) applies over the entire gate voltage range 
and confirms zv ~ 3/2. The strong gate voltage de- 
pendence of the amplitude i?o + (V g ) uncovers again the 
limitatios of one parameter scaling. It appears to be re- 
lated to the inverse of the zero temperature correlation 
length. Indeed close to the QSI transition we observe 
consistency with 



gc | 



1 1/2 



oc T, 



1/3 



(43) 



R + (V g ) - R c oc\V g — Vc 
It suggests that 

i?o+ (Vg)-R c oc 1/C(T = 0) oc \V g - V gc f oc T 1/z , (44) 

with V ~ 1/2 and z ~ 3 applies, and with it these criti- 
cal exponents apply on both sides of the QSI transition 
(see Eq. ©). 

The full quantum scaling plot R(T, V g )/Ro± (V g ) 



VS. \Vg ~ Vg C 



3/2 



/T is shown in Fig. 1101 The upper 
branch stems from the Nl-crossover and the lower one 
from the NS-transition. The two branches merge at the 
QSI-transition as the solid and dashed curves indicate. 
The solid line marks the BKT form G_ (y) given by 
Eq. ([H])) while the dashed line is G+ (y) (Eq. gTJ) 
with / = 0.065 K _1 / 3 V _3//2 and the R 0+ (V g ) estimates 
given in Fig. |9Jd. The strong gate voltage dependence 
of i?o+ exposes clearly that two parameter scaling (y, 
i?o+) is required to enter the quantum critical regime 
of the variable-range hopping controlled insulator. This 
findings are fully analogous to those observed in the 
BKT-counterpart where i?o- (V g ) exhibits a strong gate 
voltage dependence (Fig.[3fc) as well. The piecewise data 
collapse on the respective scaling function line extends 
roughly over two decades of the scaling argument. For 
large values of the scaling argument the deviations stem 
from the finite size effect with a limiting length which 
increases with reduced gate voltage (Fig. [5^,). In the 
lower branch and small scaling arguments they reveal 
that BKT fluctuations no longer dominate. The scaling 
function of the upper (Eq. I|4ip) and lower branch (Eqs. 
(123)) and (|3"Tj) ) are very different but both exhibit at the 
QSI-transition an essential singularity. 



3 Summary and Conclusions 

Over the past twenty years the analysis of sheet resis- 
tance data taken on superconducting films, expected 
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Fig. 9 a) R(V g ) /R + (V 3 ) vs. (T (V 9 ) /T) 1/3 derived from 
the sheet resistance data of Parendo et al. 6,7 . i?o+ (V 3 ) 
and To {V g ) are chosen to minimize the error in data col- 
lapse onto the scaling function exp(To {Vg) /T) 1//3 . b) Re- 
sulting estimates for -Rq+ {Vg) and To {V g ). The dashed line is 



To{V g ) ~/ 3 |V 9 -V £ 



with / = 0.065 K" 



l/3 V -3/2 and 

2.4 kQV" 1 / 2 



zv = 3/2 and the solid one i?o+ (Vg) = R, 
{V g - V 3C ) 1/2 c) R 0+ vs. T The solid line is R 0+ = R. 

,1/3 



37 kfiK" 1 / 3 ■ T f 
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-V r 2 /T(V m /K) 

Fig. 10 Quantum scaling plot R{T, V g )/Ro± (V 9 ) vs. 
\V g - V gc \ 3/2 /T. The lower branch taken from Fig. \7\ origi- 
nates from the NS-transition where Ro{V g ) = Ro-{V g ). The 
solid line marks the BKT quantum scaling function given by 
Eqs. I|2ip and (J25J). The upper branch stems from the NI- 
crossover. The dashed line is the scaling form Eq. I|40l) with 
/ = 0.065 K -1/,3 V -3 / 2 . It is characteristic for a variable- 
range hopping controlled insulating ground state and deter- 
mines the crossover temperature To {V g ) (Eq. 14111 ). 



to undergo a thickness, magnetic field, or gate volt- 
age tuned QSI transition, was predominantly based on 
the one parameter scaling form ([2(1]) . zv was usually 
determined by minimizing the error in data collapse 
using the critical sheet resistance R c and critical tun- 
ing parameter determined from the crossing point of 
the sheet resistance isotherms 4,6 . Though one pa- 
rameter scaling is asymptotically correct its applica- 
tion requires that the quantum critical regime is at- 
tained. However the generically observed saturation of 
the sheet resistance at low temperatures points to a 
finite size effect which makes this nearly impossible. In- 
deed our finite size scaling analysis of the BKT behavior 
clearly revealed that the saturation at low temperature 
is fully consistent with a limiting length L prevent- 
ing the correlation length to grow beyond L. Noting 
that disorder is irrelevant at the BKT-transition L is 
attributable to the limited homogeneity of the films. 
In agreement with previous studies [161117] we noticed 
that L decreases substantially by approaching the QSI 
transition. Thus variation of the gate voltage does not 
affect the carrier density but the inhomogeneity land- 
scape as well. Of course, improved homogeneity of the 
films would allow to circumvent these difficulties. Here 
we pursued a different route, we invoked two param- 
eter with the variable y oc \V g — V gc \ zl/ /T replaced by 
T c {V g ) jT along the BKT-transition line and T (V g ) /T 
along the insulator crossover line. The tuning param- 
eter dependent amplitude Rq± {V g ) is taken as second 
variable. At QSI criticality where Rq± (V gc ) — R c and 
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T c (V g ) oc T (V g ) oc \Vg — V gc \ zv the standard quan- 
tum scaling form is recovered. But given a strong gate 
voltage dependence of Rq± (V g ) as encountered in the 
Bi film, the outlined two parameter scaling is indispens- 
able to produce scaling plots yielding realistic estimates 
for zV. Noting that i?o- (V g ) — R c is controlled by the 
amplitude of the zero temperature correlation length 
£ (T = 0, V g ) = eo± (T = 0) \V g /V gc - 1 1 ~ v in terms of 
Bo- {Vg) - i? c oc £r7- (T = 0) (Eq. ® ) and i? + (V 9 ) - 
R c oc Co+ ( T = 0) ( E q- 0)) the strength of quantum 
fluctuations increases with reduced £o± (T = 0), a sub- 
stantial gate voltage dependence of i?o± (K,) — R c ap- 
pears to be a generic feature of the QSI-transition. Ac- 
cording to this the two parameter scaling outlined here 
opens a window to attain and enter the critical regime 
of the QSI transition and to obtain estimates for the ex- 
ponents zV, z, v and not for zv only. Indeed the value 
of V is essential to clarify the relevance of disorder at 
zero temperature via the Harris criterion [101111] and 
given z the vlidity of hyperscaling can be verified. We 
are hopeful that the two parameter scaling analysis of 
the sheet resistance will deepen our understanding of 
the QSI transition in superconducting films. 
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